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Abstract: 

It follows on general grounds that the 9 dependence in QCD is more complicated than 
suggested by the large N c approach or instanton arguments. Generically, the vacuum 
energy E vac {9) is a multi- valued function of 6 admitting the existence of metastable states. 
We discuss decays of such metastable vacua in the theory with and without the axion, and 
point out the potential relevance of this and related phenomena for constraining a dark 
matter axion. Based on the analysis of the axion potential, an idea for a new axion search 
experiment at RHIC is suggested. It is noted that the false vacuum decay proceeds with 
maximal violation of CP, even if 9 = 0. We further speculate that the famous Sakharov 
criteria for baryogenesis could be satisfied at the QCD scale. 



1 Introduction 



The axion is perhaps the most well-motivated dark matter candidate ever since the early 
1980's. Originally, the axion []J and its invisible modifications || has been introduced 
into the theory not as a dark matter candidate, but rather as the particle which solves 
an internal fundamental QCD problem -the strong CP problem. Twenty years later we 
still do not know a better solution of this problem, therefore, we should admit that the 
axion solution of the strong CP problem has successfully passed the timing test 0]. In this 
respect, unlike many other exotic particles, the axion is the unique dark matter candidate. 

In the standard thermal scenario, cosmic axions can be created through the radiative 
decay of the axion strings |4j] or due to the 'misalignment' effect at the QCD phase 
transition, see e.g. 0. These scenarios do not exhaust all possible mechanisms for the 
axion production. In particular, one more possible scenario for the axion production 
will be suggested below in the text. We note that constraints on the axion mass (or the 
coupling constant f a ~ m~ 1 ) are very sensitive to the specific scenario for production of the 
cosmic axions, inflationary model, ratio f a /T re h (where T re h is the reheating temperature 
of the universe at the end of inflation), etc. We restrict ourselves by considering possible 
constraints which follow from the QCD part of the problem. 

We start with recalling some popular ansatze for the axion potential. In the 'mis- 
alignment' mechanism, after the axion mass switches on at the QCD scale, the axion field 
begins to coherently oscillate about the minimum at 9 = 0. Such a picture corresponds 
to the simplest possible choice for the axion potential V(a) = m 2 a a 2 /2 with possible an- 
harmonic corrections, and does not take into account the periodic properties of the axion 
field a — > a + 2Trf a . One can do better and consider potential V(a) ~ cos(a// a ) motivated 
by instantons, which is a periodic function of a. One can calculate anharmonic correc- 
tions for fields near the top of the potential ||, and account for topological effects due 
to the appearance of domain walls M for this potential. All these effects are potentially 
extremely important, and can drastically change the whole picture of the cosmic axion 
production. 

However, the above potentials V(a) are the model expressions not derived from QCD. 
Moreover, they do not satisfy some general requirements of the theory, see Sect. 2, 3. There- 
fore, the corresponding calculations should be considered as qualitative estimates of pos- 
sible phenomena which may happen in the development of the axion field. 

We recall at this point that there is the one-to-one correspondence between the form of 
the axion potential V(a) and the vacuum energy E vac {9^ as a function of the fundamental 
QCD parameter 9. Indeed, the axion solution of the strong CP problem suggests that 
9 parameter in QCD is promoted to the dynamical axion field 9 — > a(x)/f a , and the 
QCD vacuum energy E vac (9) becomes the axion potential V(a). Therefore, our problem 
of analysing V(a) amounts to the study of E vac {9) in QCD without the axion. 

Recently, there has been a progress in understanding the general properties of the 
9 dependence in QCD [[?) based on a new development in supersymmetric (SUSY) 
theories (Sect. 4 and 5). The purpose of this paper is to discuss the implications of this 
new understanding for the study of local and global properties of the axion potential. As 
for the former (Sect.6.1), we note that the temperature dependence of the axion mass 
(and of entire axion potential) can be related with that of the QCD vacuum quark and 
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gluon condensates whose temperature dependence is understood (from lattice or model 
calculations). Yet, the most interesting phenomena are related to the global properties 
of the axion potential (Sect. 6. 2). As will be discussed below, the vacuum energy E vac (9) 
(and consequently, the axion potential V(a)) is generically multi- valued, i.e. in some 
parametrical regions of 9 extra local minima with higher energy could exist. It may lead 
to the phenomenon of the false vacuum decay through bubble nucleation. This process 
may provide a new mechanism for the axion production, as, when the axion is present, 
it can be supplemented by rolling to 9 = with emission of axions. Potentially, an 
account of these and other related phenomena may lead to the necessity to reconsider 
the constraints on a dark matter axion. Based on the analysis of the axion potential, we 
further suggest an idea for the new axion search experiment at RHIC (Sect. 7). Last but 
not least, we note that because CP is 100 % violated in such false vacuum decays, the 
famous Sakharov criteria || for baryogenesis could be satisfied at the QCD scale (Sect. 8). 

2 What is known about E vac (6)l 

As we already mentioned, our problem is reduced to the analysis of the 9 dependence 
of the vacuum energy due to the exact correspondence 9 —* a/ f a , E vac {9) — ► V(a/f a ). 
What do we know about E vac {9)l We know a few exact statements: 

1. We know the exact Ward Identity (WI) which (near the chiral limit m q — > with equal 
masses ) takes the form flTUfl : 

i J dx(0\T[^-GG(x) ^GG(O)} |0> = = ^(0|M|0) + O^). (1) 

Therefore, at small 9 the vacuum energy is fixed E(9) — E(0) ~ f 2 m 2 9 2 and, thus the 
axion mass m 2 a is also fixed: V(a) — V(0) ~ f 2 ml(a/ f a ) 2 , m 2 a ~ f 2 m 2 /f 2 . As it should 
be, in the chiral limit the vacuum energy does not depend on 9. 

2. The vacuum energy E vac {9) is a periodic function of 9, i.e. E vac (9 — > 2n + 9) = E vac {9). 
This periodicity is a direct consequence of quantization of the topological charge in QCD: 

|0> = 5>*"|n>, |0 + 27T> = |0>, (2) 

n 

where |n) is the winding state, n is an integer. We note that the vacuum energy in gluo- 
dynamics (QCD without quarks) is also a 2ir periodic function. 

3. In gluodynamics with the large number of colors N c the physics should depend on 9 
through the combination 9/N c in order for the U(l) problem to be solved [ |TT|] . 

4. The anomalous WI's require |12[ that Goldstone bosons which are described by the 
unitary matrix JJ^ should appear (apart from the mass term) only in the following com- 
bination with 9: 

9-iTrlogU. (3) 

This is the only allowed combination for arbitrary N c and Nf. This fact is a consequence 
of the transformation properties of the Goldstone fields U — > exp iaU and the 9 parameter 
9 — > 9 + Nf<y under the U(l) chiral rotations. 
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3 Di Vecchia-Veneziano-Witten (WW) Solution 



Based on the results outlined in the previous section, and also on the large N c argu- 



ments, Di Vecchia, Veneziano and Witten |L2| suggested the following effective potential 
to analyze the 9 dependence of E vac {9): 

W vvw {6, U) = - ^ 2 \ YM (0 - % log Det Uf - \Tr {MU + M + U + ) + ... . (4) 

In this formula M = diag(m i \(ty' l ty' l )\) and the coefficient in front of the combination 
(6 — « log Det U) 2 is the topological susceptibility in pure YM theory rather than in real 
QCD: (u 2 ) YM = < 0. 

The minimization of the potential Www (9, U) with respect to the U fields gives the 
9 dependence of the vacuum energy. The obtained result has the required 27r periodicity 
and actually is a more complicated function than the simple ansatz cos(#) ~ cos(a// Q ) 
exploited in most calculations with the axions. 

In spite of the great simplicity and attractiveness of the WW solution, this scenario 
can not be quite complete. The argument is that pure YM theory, not only QCD, should 
also obey the 2n periodicity law in 9. For YM theory in the large N c limit the potential 
(|) would imply that the 9 dependence in gluodynamics is given by Wym{9) — Wym{9 = 
0) = — {v 2 )ym9 2 /2, which can not be correct for an arbitrary 9. The first guess would 
be that our failure to reproduce the 2n periodicity is related to shortcomings of the 
1/N C expansion: the next terms, suppressed by factor 1/N C , presumably should recover 
the periodicity. However, this simple guess is not working. Indeed, any function like 
Wym(9) = — ^g-feilzM cos(9/N c ) which depends, as it should, on the combination 9/N c , 
would have period 2ttN c rather than the required 2tt. In addition, the approach of Ref. [HJ 
deals from the very beginning with the light chiral degrees of freedom and explicitly 
incorporates the Ua(1) anomaly without restriction of the topological charge to integer 
values. 

Therefore, something is missing... To understand exactly what is missing, we should 
have a better understanding of gluodynamics and its periodic properties. Only after that 
we can come back to QCD. As the first step in this direction we would like to learn some 
lessons about the 9 dependence which supersymmetric (SUSY) models offer to us. 



4 Lessons from SUSY theories 

In this section we would like to overview some SUSY models with an emphasize on the 
properties of vacuum states and their 9 dependence. As we already learned in the case 
of QCD, this is the most relevant information for the axion physics. Instrumental for 
this analysis is the 'old-fashioned' effective Lagrangian approach, in which the effective 
Lagrangian is defined as the Legendre transform of the generating functional for connected 
Green functions. Only the potential part of this Lagrangian can be fixed in this way as 
it corresponds to zero momentum n-point correlation functions. The kinetic part is not 
fixed in this framework. Thus, such an effective Lagrangian is useless for calculating the 
S'-matrix, but is perfectly suitable for addressing the vacuum properties of the theory. 
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Fortunately, this is exactly the information we need: the 9 dependence of the vacuum 
energy is the problem amenable to a study within this framework. 

To be more specific, let us consider the effective Lagrangian for supersymmetric QCD 
(SQCD). The potential part of the Lagrangian, which is fixed by the anomalous Ward 
identities, is given by [p~3|j: 



WsQCD = , 



Tr(mU) 



(5) 



where S is the gauge chiral superfield, A is the fundamental scale parameter of the theory, 
and U) = Q l Qji i,j = 1,..N* is the matter superfield with the mass matrix m. This 
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Lagrangian should describe the vacuum structure of the theory. As was found in 
there are N c different vacua in SQCD, which are labeled by the 6 angle and the discrete 
parameter k = 0, 1, ...N c — 1 such that the gluino condensate depends on these parameters 
as (AA) ~ exp ( ™t±2ik ) _ Therefore, when 9 varies continuously from to 2n, N c distinct 
and disconnected Bloch type vacua undergo a cyclic permutation: the first state becomes 
the second one, and so on. All physical quantities are periodic in 9 with periodicity 2tt, as 
these vacua can be just relabeled by the substitution k — > k — 1 after the shift 9 — > 9 + 2tt, 
keeping the physics intact [14" . 

Therefore, in SUSY theories we do have the property described in Sect. 2: the 9 de- 
pendence comes through the combination 9/N c , but the physics is 2n periodic due to 
the existence of the additional states. In SUSY models they are degenerate, in non- 
super symmetric theories we expect that they have different energies. 

Although the potential (||) has some appealing features, it can not be complete as 
was recently argued by Kovner and Shifman 0. Indeed, the scalar potential ([5]) is not a 
single-valued function of the field. If one starts at some S — So and travels continuously 
to the different (but physically equivalent) S' Q = e 2m S , the value of the potential at S' 
will be different from that at Sq. Secondly, the discrete symmetry inherent to the original 
theory is not reflected in Eq.(^). Both these shortcomings were successfully cured in || 
by the prescription of summation over all branches of the multivalued potential @. In 
the simplest case of supersymmetric gluodynamics (Nf = in Eq.(|5])), this prescription 
leads to the following definition of the effective potential Wks'- 



exp(-W KS ) =J2 ex P 



-Sin 



A N C 



+ h.c. + -(S - S)(9 + 2ixn) 
3 



(6) 



where for simplicity we suppressed the integrals over the superspace coordinates. In terms 
of the original theory, such prescription means the summation over all topological classes 



Z [ DAexp \-S + i f d i x(9 + 27rn)-^GG 

'J L J 32tt 2 



(7) 



which imposes quantization of the topological charge J d A x(l/327r 2 )GG in integer units. 
Therefore, the summation over all topological classes does not change our theory, but 
simply introduces an overall (infinite) factor into Z, which is irrelevant anyhow. It is 
clear that all anomalous Ward identities, as well as the dynamical part of the effective 
lagrangian (the first term in Eq.([])) are kept intact by this prescription. 
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It is not the purpose of this paper to discuss the technical details of this new devel- 
opment in gauge theories. We refer the reader to the original papers || [7j for details. 
Here we would like to end up this section by emphasizing some lessons to be learned from 
SUSY theories, which are relevant to us and important for the axion physics: 

1. 9 — ► 9 + 2tt is the explicit symmetry of the theory because it can be compensated for 
by a shift in n —>■ n + 1 such that the potential (^ is unchanged. 

2. The vacuum states are classified by two parameters \9, k = 0, ..N c — 1). All these states 
are physically equivalent, but they are different vacuum states labeled by the condensate 
(AA) ~ exp ( jf nk )- This additional discrete quantum number k can not be found from 
the analysis of the symmetries of the original lagrangian because of its quantum origin 
related to the anomaly. It can be understood only from explicit dynamical calculations. 

3. In SUSY models these extra states are degenerate in energy. However, in non- 
super symmetric case one could generally expect that states with different energies (i.e. 
metastable vacua) appear for some values of 9, which decay by tunneling. For the axion 
physics, this implies that decays of such metastable states may be supplemented by emis- 
sion of the axions. 

4. In general, there are domain walls in the system connecting different vacua. For the 
axion physics it would imply the existence of different kind of domain walls along with 
the ones studied previously [0 . 



5 Effective Lagrangian and 9 dependence in QCD 



In this section we describe the effective potential in QCD |7j which allows one to analyze 
the 9 dependence of the ground state. In this approach, the Goldstone fields are de- 
scribed by the unitary matrix Uij corresponding to the 75 phases of the chiral condensate: 

(*i*i) = -K* A)i^with 



U = exp 



r- 7T a A a 2 7/ 

iV2— — + i 



fn 



Nt fi 



f jrr 



where A a are the Gell-Mann matrices of SU(Nf), 7t a is the pseudoscalar octet, and f m = 
133 MeV. In terms of U and the "glueball" fields h, h, the QCD effective potential can 
be constructed analogously to the SUSY case, and has the explicit 2ti periodicity: 



-iVW(h,U) 



+00 q—1 
n=-oo k=0 



-iVW d {h, U) + i7rV I k + 



q 9 + 27m \ h — h 



P 



2tt 



2i 



(9) 



Here the "dynamical" part Wd(h, U) of the anomalous effective potential is 



W d (h,U) = ]hLog 



h \ 



p/g 



2eE 



DetU 



-TrMU +h.c. 



(10) 



where V = J d A x is the 4-volume. All dimensional parameters of this potential are 
expressed in terms of the QCD vacuum condensates, and are well known numerically: 
M = diag(mi\(^/ l ^/ l )\); the constant E is related to the QCD gluon condensate E = 
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(ba s /(32ir)G 2 ). It is interesting to note that the whole structure of Eq.(^) is rather simi- 
lar to that of the SUSY effective potential (H). Namely, it contains both the "dynamical" 
Wd and "topological" parts (the first and the second terms in the exponent, respectively). 
The "dynamical" part of the effective potential @ is similar to Wsqcd © while the 
"topological" part is akin to the improvement || of this effective potential. 

The only unknown parameters in this construction are the integers p, q, which play 
the same role as discrete integer numbers in SUSY theories, see Sect. 4. These numbers 
are related to a discrete symmetry which is a remnant of the anomaly, and can be found 
only by explicit dynamical calculations. One can argue |16| that q/p = 8 /3b where 
b = (ll/3)N c — (2/3)Nf is the first coefficient of the /^-function. However, for more 
generality we prefer not to fix q/p in what follows with the only constraint that at large 
N c -i q/p ~ l/^Vc in order for the U(l) problem to be solved. 

The heavy "glueball" fields h, h can be integrated out in Eq.@. The resulting effective 
chiral potential is periodic in 9 and takes the form |7J 



W, 



QCD{ 



U,U' 



lim — log < exp 

V->oo V \i 



VE cos 



P 



--VTr{MU + M + U + ) 



2ir 

i log Det U) + —l 
P 



/ = 0,l,...,p-l (11) 



It was argued in |7[] that Eq.(|Tl|) represents the anomalous effective Lagrangian realizing 
broken conformal and chiral symmetries of QCD. The arguments are the following: (a) 
Eq. 



LlJ) correctly reproduces the WW effective chiral lagrangian [12] in the large N c 
limit; (b) Eq.( |TTJ) reproduces the anomalous conformal and chiral Ward identities of 
QCD. 

(a) For small values of (8 — Hog Det U) < n/q, the term with / = dominates the 
infinite volume limit. Expanding the cosine (this corresponds to the expansion in q/p ~ 
1/N C ), we recover exactly the WW effective potential (|j) at lowest order in l/N c , together 
with the constant term — E = —(ba s /(32ir)G 2 ) required by the conformal anomaly: 



vvw\ 



-E- 



{v 2 )ym 



-i log Det Uf - -Tr (MU + M + U + ) + . 



'12) 



where we used the fact that at large N c E(q/p) 2 = —{u 2 )ym is the topological suscepti- 
bility in pure YM theory. Corrections in 1/N C stemming from Eq.(ll) constitute a new 
result. 

(b) It is easy to check that the anomalous chiral and conformal WI's are reproduced by 
Eq. (ITT]) . As an important example, let us calculate the topological susceptibility in QCD 
near the chiral limit. For simplicity, we consider the limit of SU(Nf) isospin symmetry 
with Nf light quarks, rrij ^qcd- For the vacuum energy for 9 < njq we obtain 



E mc (9) — —E + mJV,<0|tftf |0) cos 



(13) 



Differentiating this expression twice in 8, we reproduce the famous WI (]T|). 

We note that in general (Eq.([l3D is the particular example) the 9 dependence comes 
in a combination 9/N which naively does not provide the desired 27r periodicity for the 
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physical observables. At the same time, this formula was derived from Eq.flTTp which is 
perfectly 2tt periodic. How could it be? The answer is: the thermodynamic limit V — *■ oo 
is performed for a certain value of 9, such that only a term of lowest energy survives in 
(PI), while all other states have higher energies and therefore drop out. At the same time, 
the values 9 and 9 + 2ir are physically equivalent for the whole set of states and not for a 
selected individual vacuum state. Thus, the fact that the 9 dependence in usual V = oo 
formulas comes in the combination 9/N has nothing to do with the problem of periodicity 
in 9, as those formulas refer to one particular state out of this set [|7J. 



6 Applications to the axion physics 
6.1 Local properties of the axion potential 

Now we are ready to apply our results to the axion physics. The axion potential, by 
definition, is obtained from Eq.([0]) by the replacement 9 — > a/f a : 

W(a, U, U+) = W QCD (9 = a/ f a , U, U+). (14) 

The next step is to integrate out the Goldstone fields U exactly in the same way as was 
done in obtaining Eq . (|l"3"D for small values of 9. Technically, this problem amounts to the 
minimization of the potential W(a, U, U + ) with respect to U for a fixed value of a. 

Before discussing the general global properties of this potential, we would like to 
emphasize that all parameters of the potential (except for the integers p, q) are fixed at 
zero temperature: (M) r=0 ^ -(250 Mel/) 3 ; (a s /irG 2 ) T =o ^ 1.2 • HT 2 GeV*; f w (T = 
0) = 133 MeV. The dependence of these parameters on temperature is also (at least, 
qualitatively) known. In particular, the axion mass, which is defined as the quadratic 
coefficient in the expansion of the function E vac {9) at small 9, is proportional to the 
chiral condensate: m 2 a (T) ~ m q (0\^/^/\0)T/ f%- Therefore, m 2 (T) is known as long as 
(0|\I/\I/|0)t is known. This statement is exact up to the higher order corrections in m q . 
We neglect these higher order corrections everywhere for T < T c (T c ~ 200 MeV is the 
critical temperature), where the chiral condensate is nonzero and gives the most important 
contribution to m a . For the particular case Nf = 2 one expects the second order phase 
transition and, therefore, m\ ~ ^(0|M|0) ~ \T e -Tf for T near T c ~ 200 MeV. This is 
exactly where the axion mass does "turn on". The critical exponent in this case f3 ~ 0.38, 
see e.g. recent reviews [17] for a general discussions of the QCD phase transitions. 



In what follows we make the very plausible assumption that the temperature depen- 
dence for all observables for T <T C comes entirely through the parameters involved in the 
effective potential (the quark and gluon condensates and /„-). This assumption is based 
on the fact that the only relevant degrees of freedom at T < T c are the Goldstone bosons. 
With this assumption, the entire potential W(a, U, U + )t for T < T c (not only the first 
term of its expansion ~ m 2 ) is also known. Therefore, the procedure described above 
allows us to construct the axion potential everywhere. In particular, for small fields a in 
the limit of SU (Nf) isospin symmetry with Nf light quarks, the potential is: 



7(a) = -£7 + m ff JV>(0|**|0) T cos( T ^ r )+0(m?) > ^ < tt, T <T c . (15) 

\ Ja^f ) J a 
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Up to now we discussed only the local properties of the axion potential. Its global (or 
topological) properties may be even more important for the axion physics. 



6.2 Global properties of the axion potential 

The topological properties of the potential crucially depend on the integer numbers q,p 
we introduced earlier. Therefore, we consider the cases q — 1 and q ^ 1 separately. We 
start from the simplest q = 1 case. In this case the global properties of the potential 
qualitatively are similar to those of the WW potential (§). The only difference from 
the latter is the presence of cusp singularities at certain values of the fields, which are a 
remnant of the topological charge quantization in the effective Lagrangian approach J?|. 

An interesting qualitative phenomenon which follows from the analysis of the potential 
(|) or (PI) is a possible appearance of additional local minima at 6 ~ tt, depending on 
temperature and quark masses. Therefore, the axion potential in this region may become 
a multi-valued function, i.e. there would be two different values of the axion potential 
^1,2(0 = a /fa) for a fixed 6, which differ by the phase of the chiral field. 

In particular, for the WW potential (|j) at T = for three flavors with equal masses 
and 6 = tt, there are two degenerate states Jl2| separated by the domain wall. The wall 
surface tension in this case was recently calculated by Smilga fllB]: 



o- = 3V2 M - j m«fl (16) 

For 9 7^ 7r the energies of the vacua are not degenerate anymore but are splitted apart 
by the amount AE ~ m q (6 — it). For the case of equal masses, a metastable state in the 
WW scenario exists in the region tt/2 < 9 < 3tt/2 fl8| . (As we will discuss shortly, this 
phenomenon becomes far more general for the potential (ffl|) if q ^ 1.) 

A metastable vacuum decays to the ground state with the formation of bubbles of the 
stable phase. The quasiclassical formula for the decay rate per unit time per unit volume 
was derived many years agoQ [|T9| : 

r oc exp (- ^E) 3 ) " ex P(-^4)- (17) 
Using (0), one finds that in the WW scenario the lifetime of the metastable state at zero 



temperature is much larger than the age of the Universe ||18j| . Such phenomenon might 
play an important role in the development of the early Universe during the QCD epoch. 
However, to make the corresponding estimates one can not literally use Eq.(|D]) because 
the parameters which enter this formula depend on temperature and may drastically 
change the result near the phase transition point. It is quite possible that during the phase 
transition the relevant factor which enters ([l?]) vanishes: a A /(AE) 3 ~ \T C — T\ a — > 0, a > 
at \T C — T\ —>■ and, therefore, T ~ 1. Indeed, in the mean field approximation where 
(0|M|0) T ~ (T c - T) 1//2 , fl ~ (T c - T) we have a ~ (T c - T) 3 / 4 and AE ~ (T c - T) 1 ' 2 . 
Therefore a 4 /(A£) 3 ~ \T e - T| 3 / 2 -> 0. 

1 This formula is derived for zero temperature, for a very high temperature [ p0[ the correct expression 
is T oc cxp{-167rcr 3 /(3T(Ai;) 2 )}. 
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Now we assume that the axion exists. In this case, if the axion field at temperature 
T ~ T c is trapped in a metastable minimum at 9 ~ tt, there are two different options for 
the decay of this false vacuum. If T ~ 1 at T ~ T c , tunneling is not suppressed, and the 
false vacuum tends to lower its energy first by bubble nucleation without changing 9, and 
then by relaxing to 9 = with production of axions. This is advantageous since the axion 
potential is nearly flat. However, if T is still very small at T ~ T c , the main mechanism 
for the decay could be related to the direct axion production. 

For the case q ^ 1 (which we prefer, see ||16||), the physics is even more interesting. 
Additional metastable vacua appear for arbitrary quark masses and for an arbitrary 9, 
not necessarily in the region 9 ~ tt. Therefore, decays of these false vacuum states will 
be a general phenomenon. In this case the wall surface tension can be easily calculated 
from ( |TTD and for 9 = is given by (see [21] for details) 




+ 0(mX 



This formula substitutes Eq.flTB"!) describing the wall surface tension at 9 ~ tt for the 
WW potential (|). A distinct difference between these two cases is the absence of the 
chiral suppression ~ m q in Eq.(^), which apparently would make penetration through 
the barrier even more difficult in comparison to the WW potential. The energy splitting 
between the ground state and metastable state at 9 = is 



AE = m q N } 



2tt 



1-cos — - +0(mf), (19) 



qN f 



and we obtain 

3 3 • 2 7 • ttV f%E* (l ~ cos |) 4 27 ir*<?N f f^G 2 f 



S, 



q±Nf M3 fi_ cos ^ 3 256 p* m 3 



(20) 



qN f J 

Eq.(pO|) shows that the parametric suppression of the decay is largely overcome due to 
a numerical enhancement. The latter depends crucially on the particular values of the 
integers p,q. In particular, for our favorite choice p = 11N C — 2Nf,q = 8 fl6fl , Eq. (|20|) 
yields a factor ~ 10, while e.g. for p = N c ,q = 1 (as motivated by SUSY, see pif ) it is 
approximately two orders of magnitude larger, but still much smaller than the estimate 
of |TS|| for the WW potential. Of course, all remarks made above concerning the temper- 



ature dependence of such effect apply for the estimate (|20|) as well. In particular, T could 
be of order one during the QCD phase transition in development of the early Universe 
and, therefore, a metastable state could decay just before the Universe cools down to the 
temperature where Eq. (|20|) could be applied. 

Another related question is the dynamics of the domain walls which separate two 
different vacua with the small splitting AE ( |19|) . These domain walls are absolutely 
harmless for development of our universe because they decay in proper time scale. Indeed, 
the required value for the pressure (which is equal to AE in our notations) for the safe 
decay of the wall is given by |22j 

2tt \ ^ 
qNf) ~ MV 



AE = m q N f 



(W) ( 1-cos — ] (21) 
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where Mp is the Plank scale. Inequality Q2T| ) is perfectly satisfied with our parameters 



AE and a Therefore, such domain walls do not lead to a cosmological disaster, 



but rather may have very interesting cosmological consequences which still have to be 
explored. Similar domain walls which separate vacua with different phases of the gluino 
condensate have been recently discussed by Shifman in SUSY models |23|. It is quite 



remarkable that this phenomenon may exist not only in SUSY models but also in the 
physically relevant case of QCD. 

We conclude this section by emphasizing that the phenomena just described may lead 
to the necessity to reconsider the constraints on a dark matter axion because the dynamics 
of the axion field could be more complicated than it was originally thought. We stress 
that the decay of the metastable state described above proceeds (with or without axions) 
through the expansion of bubbles with 100% violation of CP invariance. This is because 
the phase of the chiral condensate in the metastable vacuum is nonzero and of order 1. 
This leads to violation of CP even if 9 = 0. (This is not at variance with the Vafa- 
Witten theorem [24|| which refers to the lowest energy state only.) It may have profound 
consequences for the development of the early Universe at the QCD scale because such 
effects could lead to a new mechanism for baryogenesis! Indeed, the famous Sakharov 
criteria |J could be satisfied in this scenario, see Sect. 8 for a more detailed discussion. 



7 New axion search experiment at RHIC? 

The development of the early Universe is a remarkable laboratory for the study of most 
nontrivial properties of the particle physics. What is more amazing is the fact that 
these phenomena at the QCD scale can be, in principle, experimentally tested at RHIC, 
Brookhaven. 

We expect that, in general, an arbitrary |#)-state would be created in the heavy ion 
collisions, similarly to the creation of the disoriented chiral condensate (DCC) with an 
arbitrary isospin direction. It should be a large domain with a wrong 9^0 orientation. 
As we shall see in a moment, for both cases (DCC and |#)-state) the difference in energy 
between a created state and the lowest energy state is proportional to m q and negligible 
at high temperature. Therefore, energetically an arbitrary \9) can be formed. The way 
of how the created |#)-state will relax to the ground state of lowest energy \9 = 0) is a 
separate issue. If, somehow, an equilibrium state with a large correlation length in the 
large volume V is formed, the only possible way to relax 9 to zero would be the axion 
mechanism (if they exist). However, due to the fact that we do not expect to create 
an equilibrium state with an infinite correlation length in the heavy ion collisions, the 
decay of a |0)-state will also occur due to the Goldstone U fields with specific CP-odd 
correlations^. Therefore, two mechanisms of the relaxation of a |#)-state to the vacuum 
would compete: the axion one and the standard decay to the Goldstone bosons. In the 
large volume limit if a reasonably good equilibrium state with a large correlation length 
is created, the axion mechanism would win; otherwise, the Goldstone mechanism would 
win. In any case, the result of the decay of a |6*)-state would be very different depending 
on the presence or absence of the axion field in Nature. 

2 A similar phenomenon has been recently discussed in Ref. pa 
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Before going into details, we would like to recall some general properties of the DCC 
which (hopefully) can be produced at RHIC (see e.g. for a review), with an emphasize 
on the analogy between the DCC and a misaligned |#)-state. If the cooling process is 
very rapid and, therefore, the system is out of equilibrium, there will be a large size of 
the correlated region in which the vacuum condensate orientation mismatches its zero 
temperature value. The absolute value of the chiral condensate right after the phase 
transition is expected to be close to its final (zero temperature) magnitude. However the 
vacuum direction of the formed condensate is still misaligned since it takes a longer time 
for the vacuum orientation to relax due to the small free energy difference ~ m q between 
the formed and true vacuum states. 

To be more specific, let us consider the case Nf = 2. The matrix U is parametrized 
by the misalignment angle and the unit vector n in the isospin space: 

U = , Tr{r a T b ) = 25 ab , UU + = 1 , (¥ L ¥ R ) = -\(% L * R )\ (22) 

The energy density of the DCC is determined by the mass term: 

= —-Tr(MU + M + U + ) = -2m\ (M) | cos(0) (23) 

where we put m u = = m for simplicity. Eq.(p3|) implies that any ^ 0(mod 2ir) 
is not a stable vacuum state because -7^|^ ^ 0, i.e. the vacuum is misaligned. On 
the other hand, the energy difference between the misaligned state and true vacuum with 
= is small and proportional to m q . Therefore, the probability to create a state with 
an arbitrary at high temperature T ~ T c is proportional to exp[— V(E$ — Eq)/T] and 
depends on only very weakly, i.e. is a quasi-flat direction. Right after the phase 
transition when becomes nonzero, the pion field begins to roll toward = 0, and 

of course overshoots = 0. Thereafter, oscillates. One should expect the coherent 
oscillations of the ir meson field which would correspond to a zero-momentum condensate 
of pions. Eventually these classical oscillations produce real 7r mesons which hopefully 
can be observed. In a sense this picture is very similar to the standard "misalignment" 
mechanism || |3|] for production of the cosmic axions during the QCD epoch. 

Now we turn to our main point when the U(1)a phase of the disoriented chiral con- 
densate is also nonzero and, therefore, the |0) -vacuum state could be formed. To take 
into account this U(1)a phase we choose the matrix Uij in the form U = diag (e 1 ^). The 
energy density of the misaligned vacuum is determined in this case by Eq. flTT|) . The most 
important difference between Eqs. (|23"D and (|TT|) is the presence of the parametrically large 
term ~ E 3> m q \(ty^) \ in the expression for energy (|TTD, describing the U(1)a phase of 
the disoriented chiral condensate. This term, as was explained in Sect. 5, is related to the 
anomalous WI's, and does not go away in the chiral limit. 

The key point is the following. For arbitrary phases 0j the energy of a misaligned 
state differs by a huge amount ~ E from the vacuum energy. Therefore, apparently there 
are no quasi-flat directions along 0j coordinates, which would lead to the long wavelength 
oscillations with production of a large size domain. However, when the relevant combi- 
nation (52i4>i — @) from Eq.flnp is close by an amount ~ 0(m q ) to its vacuum value, a 
Boltzmann suppression due to the term ~ E is absent, and an arbitrary misaligned In- 
state can be formed. In this case for any 9 the difference in energy between the true \9) 



11 



vacuum and a misaligned \9)- state (when the fa fields are not yet in their final positions 
fa(9)) is proportional to m q and very small in close analogy to the DCC case. 

After this point we can apply the same philosophy as for DCC. The chiral fields^ fa 
begin to roll toward the true solution fa(9) and of course overshoots it. The situation is 
very similar to what was described for the DCC with the only difference that in general 
we expect an arbitrary |#)-disoriented state to be created in heavy ion collisions, not 
necessarily the \9 = 0) state. The difference in energy between these states is propor- 
tional to m q , as follows from the fact that the 9 dependence of any physical observable is 
proportional to m q . 

If a reasonably stable | (9 ^ 0) state is created, it could decay into the axionsQ. To 
estimate the effect, we consider an ideal case when all energy stored in the \9 ^ 0) state 
will be released through production of the axions^. In this case from Eq.(|13D we expect 
the following axion density n a for a random parameter 9 ~ 1: 

„„ = ,rVy/|(f n U - cos ') ~ , . 10" MeV 3 ( , (24) 

where parameter p is a suppression factor, see footnote (5). The standard way to detect the 
produced axions is to use their property of conversion to photons in an external magnetic 
field [p8fl . The relevant axion-photon coupling constant is defined in the following way: 



L = g J^L aF ~ F g ( DFSZ ) = -?- g ail (KSVZ) = -^- h - (25) 

where we specified the coupling constants in the limit m u = for two popular models 
0. The conversion probability is given by p8fl : 

P*~r = \&„ i BH2 ) ( 26 ) 

where B is the external magnetic field and I its length. With our estimate (|24| ) for the 
axion density, one can get the following formula for the probability of detection of a 
photon from the axion conversion [ 27 1 : 

P 7 ~ P a ^n a V ~ pm a (—) 2 (B 2 l 2 ) V, (27) 

2,71 v ' 

where we took into account the relation /^m^ ~ mq^l^^lO), see Sect. 6.1. In this formula 
V ~ A 3 where A is the correlation length for the misaligned |#)-state. As was discussed 

3 If 9 0, the Goldstone fields are not exactly the pseudoscalar fields, but rather are mixed with the 
scalars; the mixing angle between the singlet and octet combinations also depends on 9, see |2l[] for detail. 

4 The possibility of a production of the axions in the heavy ion collisions was independently discussed 
by Melissinos |^| . AZ thanks Adrian Melissinos for a conversation on the subject. 

5 We are, of course, aware that this ideal case can not be realized in Nature. There should be a 
strong suppression factor p in front of Eq.(^) due to the fact that there is no superselection rule in a 
finite volume noncquilibrium state. Therefore, the 9 parameter can be relaxed to zero by others means, 
not necessarily related to the axion productions. The parameter p should depend, first of all, on the 
correlation length A of the formed misaligned \9) state. We do not know at the moment how to estimate 
this parameter p{\). For the ideal case, p = 1. 
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above, the mechanism of production of axions by the coherent processes corresponds 
to a condensate of the "almost" zero-momentum (~ m q /X) axions. If the collision is 
asymmetric (as was suggested in P7[| ), the photon signal from the axion-photon conversion 
will be unambiguous because it will be the "almost" monochromatic wave with Au ~ 
m q /X and u to be determined by the asymmetry in energies of the beams. 

Therefore, the heavy ion collisions give us a unique chance for a new axion search 
experiment. We would like to stop here with this optimistic note. 



8 Conclusion 

• Probably the most important outcome of our study is a better understanding of the 
global and local properties of the axion potential V(a), which may result in the necessity 
to reconsider the constraints on a dark matter axion. The temperature dependence of the 
potential at T < T c is expressed in terms of the T-dependent vacuum condensates which 
have been extensively studied on the lattice and in the models |l7j and are fairly known. 
Therefore, we know the temperature dependence of the potential as long as we know the 
temperature dependence of the chiral and gluon condensates. 

• The axion potential is generically a multi-valued function. Such a property may 
lead to the phenomenon of the false vacuum decay through bubble nucleation, which may 
make the dynamics of the axion field far more complicated (and interesting) than was 
previously thought. 

• Based on the analysis of the potential V(a) expressed in terms of the condensates 
and the knowledge of their temperature dependence, we suggested a new idea for the 
axion search experiment at RHIC. 

• As a byproduct of our analysis (which was mentioned at the end of Sect. 6), we 
would like to speculate on the possibility of baryogenesis at the QCD scale (with or 
without axions). Indeed, it appears that all three famous Sakharov criteria M could be 
satisfied in the decay of a metastable state discussed above: 

1. Such a metastable state is clearly out of equilibrium; 

2. CP violation is unsuppressed and proportional to m u mdm s 9 e ff, 9 e ff ~ 1. As is known, 
this is the most difficult part to satisfy in the scenario of baryogenesis at the electroweak 
scale within the standard model for CP violation; 

3. The third Sakharov criterion is violation of the baryon (B) number. Of course, the 
corresponding U(l) is an exact global symmetry of QCD. However, a "spontaneous" 
breaking of the baryon U(l) symmetry (as a result of interactions with the domain wall) is 
not forbidden, and would be sufficient. The latter may arise in a way similar to electroweak 



baryogenesis, see e.g. |29| for a review. 

The viability and details of such a mechanism for baryogenesis are still to be explored 
fjQfl , however, in general, one could expect a large asymmetry according to the unsup- 
pressed CP violation ~ m u mdm s / A 3 in bubble nucleation. What is amazing is the fact 
that each step in such a scenario for baryogenesis could be, in principle, experimentally 
tested at RHIC, Brookhaven. 

• The remark that QCD scale could be an interesting place to look is also motivated 



by recent observations of the MACHO collaboration [p]J . In particular, the domain walls 
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with QCD scale which always separate the true and false vacuum states do not lead to 
a cosmological disaster, see (|21|). They rather may be very interesting objects for the 
problem of the structure formation at small scales ~ M & . 
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